AD-A089  846 


UNCLASSIFIED 


BROWN  UN IV  PROVIDENCE  RI  LEFSCHET2  CENTER  FOR  DYNANI— ETC  F/8  12/1 
SPLINE  APPROXIMATION  FOR  AUTONOMOUS  NONLINEAR  FUNCTIONAL  DIFFER— ETC  ( U ) 
JUN  80  F  KAPPEL  D*A«29-79-C-«161 

AFA«U»-Tn-An-0MA9  M 


END 

IM90 


OTIC 


A0898 


UMcnTATION  page 


ffiOSR-ftl-  80.086 


r^r  messifijee —  READ  INS  t*RUCT  ONs" 

_ nF.FORE  COMPLETING  FORM 

GOVT  ACCESSION  NO.  S.  RECIPlEN T‘S  CAT  ALOG  NUMBER  ~ 


«.  TITLE  (and  Subttll9) 

SPLINE  APPROXIMATION  FOR  AUTONOMOUS 
NONLINEAR  FUNCTIONAL  DIFFERENTIAL 
EQUATIONS 

7-  AUTHORS.) 

FRANZ  KAPPEL 


5.  TYPE  OF  REPORT  »  PERI 00  COVEREO 


Interim _ 

,  S.  PERFORMING  ORG.  REPORT  NUMBER 


».  CONTRACT  OR  GRANT  NUMBERS) 


AFOSR  76-3092  ✓ 


9.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS  '0-  PROGRAM  ELEMEIJT  PROJECT.  TASK 

AIR  FORCE  OFFICE  OF  /  area  work  un  numbers 

DIVISION  OF  APPLIED  MATHEMATICS  / 

BROWN  UNIVERSITY,  PROVIDENCE,  R.I.  02912  61102F  2304/A1 

II.  controlling  office  name  AND  AOORESS 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARC 
BOLLING  AIR  FORCE  BASE  /WM  '*•  number  of  pages 

WASHINGTON.  D.C. _ [  .21 - 

TT  MONI TORIN AGENCY  NAME  6  ADDRESS^//  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (ol  thl e  report) 

UNCLASSIFIED 


IS*.  DECLASS1  MCA  T:  ON/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetrmct  entered  In  Block  20,  II  dlflerent  from  Report) 


19.  KEY  WOROS  (Continue  on  eeveree  elde  II  neceeeory  end  Identity  by  block  number) 


20.  A BST RA Continue  on  rover ee  elde  II  neceeeery  end  Identity  by  block  number) 

Based  on  abstract  approximation  results  in  semigroup 
theory,  wo  devolep  an  approximation  scheme  for  nonlinear  autonomov 
funct  lonal -dif^rejtfial  equations  with  globally  Lipschitzean  rigl 
hand  sidefrC5 Tne  sdrreme  can  be  realized  by  using  spline  approxi¬ 
mation  of  The  state. 

V 


DD  ,Fj2SMt,  1473  COITION  OF  I  NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 


AF0SRI.TR-  8  0-0863 

)  - 

•:Ji)  i 


leva-  ® 


SPLINE  APPROXIMATION  FOR  AUTONOMOUS  \ 

y  *  *  \ 

NONLINEAR  FUNCTIONAL  JDIFFERENTIAL  EQUATIONS.  ^ 


/  ^  %/Kappel^ 

1  — “"MsTCUemaFiscfies  Institut 

UniversitHt  Graz 
Elisabethstrasse  11 
A  8010  Graz  (Austria) 


UJ 

;  J 


Lefschetz  Center  for  Dynamical  Systems 
Division  of  Applied  Mathematics 
Brown  University 
Providence,  R.  I.  02912 


jJ]a 


/  * 

•cjw 

/ _ 


y  ,  , 

<4  > 


This  research  was  support e^jji  part- .-by— J^Lr  Force  Office  of  Scientific 
Research  under  contract  <^PAFOSR<»76-3092^t  <n  part-  by  the  tin-trod 
States  Army  Research  Off ice~uh3er  contract  #AROADAAG29-79-C-0161^ 


yo/f  i  7 

_  Approve’,  fr»*  i-ibllo  release, 

8  0 


SPLINE  APPROXIMATION  FOR  AUTONOMOUS 


NONLINEAR  FUNCTIONAL  DIFFERENTIAL  EQUATIONS 

by 

F.  Rappel 


ABSTRACT : 


Based  on  abstract  approximation  results  in  semigroup  theory  we 
develop  an  approximation  scheme  for  nonlinear  autonomous  functional- 
differential  equations  with  globally  Lipschitzean  right-hand  side.  The 
scheme  can  be  realized  by  using  spline  approximation  of  the  state. 
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.  Introduction  and  Notation 


In  this  paper  we  show  that  approximation  techniques  developed  in  [5] 
for  linear  autonomous  functional-differential  equations  of  retarded  type  are 
also  applicable  to  a  rather  broad  class  of  nonlinear  equations.  These 
techniques  already  have  been  successfully  applied  to  linear  autonomous 
equations  of  neutral  type  (see  for  instance  [9],  [10]).  As  in  the  linear 
case  the  approach  here  is  based  on  abstract  approximation  results  in  semi¬ 
group  theory  and  provides  a  sequence  of  ordinary  differential  systems  of 
increasing  dimension  whose  solutions  approximate  those  of  the  original  delay 
system.  The  same  algorithm  is  obtained  in  [2]  by  a  different  method  under 
the  more  restrictive  condition  of  differentiability  of  the  right-hand  side. 

On  the  other  hand  the  proof  in  [2]  also  works  in  the  time-dependent  case. 

The  general  idea  of  using  abstract  approximation  results  of  semigroup  theory 
in  order  to  get  algorithms  for  the  numerical  solution  of  optimal  control 
problems  involving  delay  systems  goes  back  to  [3],  [4].  Here  we  restrict 
ourselves  to  the  approximation  of  nonlinear  autonomous  delay  equations.  With 
respect  to  control  problems  for  nonlinear  delay  systems  see  [1].  The  assump¬ 
tions  of  [1]  in  the  autonomous  case  are  stronger  than  those  imposed  here. 

In  [11]  the  scheme  of  averaging  projections  developed  in  [4]  was  shown  to  be 
convergent  also  in  the  case  of  nonlinear  autonomous  equations  with  locally 
Lipschitzean  right-hand  side.  For  nonlinear  time -dependent  equations  satisfy¬ 
ing  Caratheodory  type  conditions  an  approximation  scheme  based  on  interpolation 
by  splines  of  first  order  was  developed  in  [12]. 

Our  state  space  will  be  Z  *  Hn  x  l^-r.O-IR11) ,  r  >  0,  which  is  a  Hilbert 
2  2  2 

space  with  norm  |  (n » 4>)  |  ^  “  I  n  |  +  1 1 2  an^  corresponding  inner  product. 


2. 


| • |  and  | * | ^  denote  the  Euclidean  norm  in  lRn  and  the  usual  norm  in 

L„(-r,0jlRn) ,  respectively.  It  will  be  necessary  to  endow  Z  with  an 
(  2  2  2 

equivalent  norm  |(n><t>)|  “In!  +  |<J>(s)|  g(s)ds,  where  g  is  a  positive 

8  J 

weighting  function  on  [-r,0].  Z  supplied  with  this  equivalent  norm  and  the 

T  x 

corresponding  inner  product  (n2><J>2^g  *  nln2  +  j  $2  (s)g(s)ds 

will  be  denoted  by  Z^.  i^(-r ,0 ;JRn)  is  the  linear  space  of  square-integrable 
functions  [-r,0]  -*]Rn  (in  contrast  to  L^C-r.OjlR*1)  which  is  the  Banach 
space  of  equivalence  classes  of  such  functions).  W  *  denotes  the  space  of 
absolutely  continuous  functions  <J>:  [-r,0]  ■+■  ]Rn  such  that  <{>  is  square- 
integrable.  5^,  k  =  0,1,...,  is  the  subspace  {  (<}>  (0)  ,4>)  1 4>  6  C^J  of  Z, 
where  C  is  the  space  of  k-times  continuously  differentiable  functions 
<J>:  [-r,0]  -*•  ]Rn.  Finally,  given  a  function  x:  [-r,a)  ->-]Rn,  a  >  0,  the 
function  xt,  t  €  [0,a),  is  defined  by  xt(s)  *  x(t+s),  s  e  [-r,0]. 

2.  The  Class  of  Autonomous  Nonlinear  Delay  Equations 

Let  a  map  fj  is!^(-r,0jRn)  -*  1R°  be  given  and  assume  that  the  following 
two  conditions  hold  throughout  the  paper: 

(HI)  If,  for  some  a  >  0,  x  is  in  L2(-r,a;IRn)  and  has  a  representation 
which  is  continuous  on  (0,a),  then  for  any  representation  )(  of  x 
the  map  defined  by  t  -*■  f (xt )  on  [0,a]  is  in  the  same  equivalence 
class  of  L^(0,a^n). 


(H2)  There  exist  constants  L  >  0  and  r  J  •  0,...,m,  with 

0  -  r_  <  r  <  ...  <  r  •  r  such  that  for  any  in  i*?(-r ,  0  jIRn) 

u  J.  in  l 

I f (♦)  -  f(*)|  <  L(  f  | <K-r ,)  -  »(-r.)|  +  |<MH2>- 
-  j-0  J  J 
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The  Cauchy  problem  we  are  concerned  with  in  this  section  is 

x(t)  «  f(xt)  ,  t  _>  0, 

x(0)  “  rj  €]Rn,x(s)  *  <J>(s)  a.e.  on  [-r,0]. 


(2.1) 

(2.2) 


where  <p  €  L£(-r  ,0*]lln) .  By  a  solution  of  (2.1),  (2.2)  we  mean  a  function 
x:  [-r,ot)  -*]Rn,  a  >  0,  such  that  (2.2)  holds  and  x(t)  =  n  +  f(x  )ds  for 

J0  s 

t  €  [0,a).  Condition  (HI)  just  assures  integrability  of  f(x  )  for  all 

s 

functions  x  which  are  candidates  for  a  solution. 


Lemma  2.1. 

a)  For  any  (q,(J>)  €  Z  problem  (2.1),  (2.2)  has  a  unique  solution 
x(t)  =  x(t;n,<)>)  existing  on  (-r,®).  Moreover,  this  solution 
is  continuously  dependent  on  initial  data. 

b)  The  family  T(')  of  operators  Z  Z  defined  by 

T(t)(n,4>)  -  (x(t;n,<j>),  xt(n,<f>)),  t  >  0,  (n,4>)  e  Z,  has  the 

following  properties: 

(i)  T(0)  =  I. 

(ii)  T(t+s)  -  T(t)T(s) ,  t,s  >_  0. 

(iii)  For  any  (n,4>)  €  Z  the  map  defined  by  t  -*■  T(t)(n.<l>) 
is  continuous  on  [0,°°). 

(iv)  There  exist  constants  M  ^  1  and  u)  €  1R  such  that 

| T(t ) -  T(t)Zj  |  <  Me^*  |  z^~z2^  f°r  0  and  zj  ^  z» 

z  z 

i  -  1*2 , 


4. 


This  lemma  is  a  special  case  of  results  proved  in  [11]  (Proposition  1.1, 
Remark  1.7,  Proposition  1.2  and  Proposition  2.1).  Note,  that  (HI)  and  (H2) 
imply  that  the  global  Borisovifi-Turbabin  conditions  as  formulated  in  [11] 
hold. 

In  [11;  Proposition  2.2]  it  was  shown  that  the  infinitesimal  generator 
Sztf  of  the  semigroup  T(*)  defined  in  Lemma  2.1  is  given  by 


dom_Q^ **  ((4>(0)  ,<J>)  |4>  e  W1’2}, 

^(4>(0),<(>)  -  <f  (♦),*)  for  WO)  A)  Edom  q/. 


(2.3) 


We  may  decompose  a/, 

&  «  _0f£)  +  (2.4) 

where  J^(<KO),<|>)  -  (0,<j>)  and  ja^(<(.(0)  ,<j>)  -  (f(<|>),0).  -Q/q  is  a  closed 

linear  operator  with  dora  ja^  *  dom  _Qf  dense  in  Z.  is  the  infinitesimal 

generator  of  the  C^-semigroup  corresponding  to  the  equation  x(t)  *  0  (cf. 

till). 

If  &  is  a  single  valued  nonlinear  operator  in  a  Hilbert  space  X, 
then  -  wl,  o)  €  1R,  is  called  dii  iipative  if 

2 

£  wlxl-x2^X 

for  all  xA  €  dom  #,  i  «  1,2.  In  general,  _Q^  -  0)1  will  not  be  dissipative 
for  any  o)  €  1R  (ja^  the  operator  defined  in  (2.3)),  because  the  constant  M 
appearing  in  (iv)  of  Lemma  2.1,  b)  cannot  be  chosen  as  one.  But  following  an 
idea  given  by  Webb  in  [15]  it  is  possible  to  renorm  Z  such  that  we  get 


5. 


dissipativeness.  Corresponding  to  the  numbers  r^  in  (H2)  we  choose  a 
weighting  function  g  satisfying 

g(s)  =  m  -  j  +  1  on  (-r^.-r^),  j  = 

Of  course,  |*|  is  a  norm  on  Z  equivalent  to  |*|„. 

s 


Lemma  2.2. 

a)  There  exists  an  u>  £  1R  such  that  srf  -  0)1  is  dissipative  in  Z^. 

b)  There  exists  a  X^  >  0  such  that 

&(I  -  \jl/)  =  Z 

for  X  £  (0,XQ). 

c)  T(t)z  =  lim  (I  -  -j^O  "z  for  all  z  €  Z  uniformly  with  respect 

n-*»  n 

to  t  in  bounded  intervals. 


Proof.  For  (4>(0)  ,4>)  and  Op(0),iJj)  in  domjaf  we  obtain  (using  (2.3)  and 
(H2)) 


<_c/(«t>(0)  ,<{))  -jtf(*( o)  ,</>),  (<|»(0)-^(0)  > 

®  — r 

m  I  I  I  I  1  5  f  J“i  t 

<  L(  l  <K-r.)-«M-r.)  +  <M> L)  <K0)h|»(0)  +  l  (m-j+1)  (*-  iP)  (d>-i|>)ds 

j*0  3  3  1  1  1  1  j-1  J-rj 


2  »  9 

2  .  m 

1  2L 

(<K0)-iK0),<M0 

+  |l2 

g 

4>(0)h|>(0) 

+  \  l  |<K-r,)-iH-r  ) 
j-1  3  3 

2  .  m 

♦  (O)hKO) 

-1  l 

j-i 

<J>(-rj)-^(-rj) 

1 

2 


2 


<  12L  +  j(L2+1)] 


(<t>(0)-^(0)  ,4>-4») 


2 

g 
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r\ 

This  proves  a)  with  u  =  2L  +  ^(L  +1) .  In  order  to  prove  b)  we  show  that 

(I  -  \30(<KO),4>)  =  (n»  (2.5) 

has  a  solution  (4>(0),4>)  £  domJ^  for  any  (n.'JO  £  Z  provided  X  is 
sufficiently  small.  (2.5)  is  equivalent  to 

<t>(0)  -  Xf  (4>)  =  p  and  <f>  -  X(f>  =  (J). 

The  second  equation  implies  (s)  =  eS^<j>(0)  ~  X  (  e^S-T^/,^(T)dT  which  for 

12  ® 

any  (|)(0)  is  certainly  in  W  *  .  For  (f)(0)  we  get  the  fixed  point  equation 
(f)(0)  =  h ((f)(0))  in  Kn,  where 

h(a)  =  Xf  (e  ^  j  *~T^\i(T)dT)  +  r). 

An  easy  calculation  using  (H2)  shows  that  h  has  the  Lipschitz  constant 
XL(m+l+r  )  on  ]Rn.  This  proves  b)  with  XQ  =  l/LCnri-l+r1^2) . 

In  Ill;  Proposition  2.3]  it  is  shown  that  u(t)  =  T(t)z  for  z  £  dom^ 
is  a  strong  solution  of  the  abstract  Cauchy  problem 

u  =  ja/u,  u(0)  =  z. 

On  the  other  hand  any  strong  solution  of  this  Cauchy  problem  is  represented 

by  u(t)  =  T(t)z,  where  T(t)z  =  lim(I  -  (see  [8;  Theorem  II)). 

n-x»  ” 

Therefore  the  semigroups  T(.)  and  T(*)  coincide  on  dom  a/.  Since  dom 
is  dense  in  Z,  we  have  T(t)  -  T(t),  t  >  0. 

We  express  the  situation  of  Lemma  2.2,  c)  by  saying  that  the  semigroup 
T(-)  is  generated  by  a/. 
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For  our  approach  we  shall  need 

Lemma  2.3.  The  sets  j^k,  k  =  0,1,...,  and  (I  -  k  =  1,2,..., 

are  dense  in  Z  for  A  sufficiently  small.  Moreover,  5^k  c:  dom  sxl  for 
k  =  1,2, . . .  . 

Proof.  Define  5*k  =  {(<fr(0),$)|$  €  Ck  and  <J>(0)  =  f (<)>)},  k  =  1,2 . 

Then  (I  -  Ao^)Ck  3  (I  -  =  S^k~^  m  dom  _Q/k  and  everything  follows 

from  the  fact  that  domJ3^k  is  dense  in  Z.  Note,  that  is  the 

infinitesimal  generator  of  a  (^-semigroup  of  bounded  linear  operators.  It 
only  remains  to  show  (I  -  Aja^)<^k  =  5^k  It  is  already  established  that 
equation  (2.5)  has  a  unique  solution  <<f>(0) ,4>)  for  (n > =  (i|»(0)  ,iJj)  , 
ip  €  Ck  1  and  A  €  (O.A^).  ip  €  Ck_1  immediately  implies  <{>  €  Ck.  From 
<t>(0)  -  A<j)(0)  =  4K0)  and  <f>(0)  -  Af(\p)  =  ^(0)  we  get  $(0)  =  f (^)  which 
proves  ((f)(0), <j>)  €  ^k.  ^k  «=  domJ3^  is  clear  by  (2.3). 

Remark.  From  Lemma  2.2,  p)  and  b)  it  follows  that  (I  -  A_£*0  ^  exists  for 
A  €  (0,Aq)  and  is  a  globally  Lipschitzean  operator  Z  dom  Ji/.  Therefore 
^k  =  (I  -  A_q^)  '*'^k  ^  for  A  €  (0,Aq),  which  implies  that  ^k  is  also 
dense  in  Z,  k  =  1,2,...  . 

3.  The  Approximation  Scheme 

Following  the  idea  already  used  in  [5]  we  choose  a  sequence  Z^  of 
finite  dimensional  subspaces  of  Z  such  that 


■  -  -  -  -  — . 
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ZN  c  dom  _Q<  N  =  1,2 . 

Let  PN:  Zg  -*■  Z^j  be  the  orthogonal  projection  onto  Z^  and  define  _Q^  by 

"^N  -  PN"^^>N  *  N  =  1 , 2 . 

We  call  the  sequence  {ZN,PN,ja^}  an  aPProximation  scheme  for  equation  (2.1). 

Lemma  3.1. 

a)  J3^  -  wl  is  dissipative  in  Z  for  all  N,  where  o)  is  the 
constant  of  Lemma  2.2,  a). 

b)  is  globally  Lipschitzean  on  Z,  N  =  1,2,,..  . 

c)  For  all  A  >  0  such  that  -  A-Q^)  =  Z 

(I  -  AJ*^)-1  ZN  c  ZN,  N  =  1,2,...  . 

Proof .  Part  a)  is  proved  in  the  same  way  as  in  the  linear  case  (see  [5] , 

Proof  of  Theorem  3.1).  For  the  proof  of  part  b)  we  recall  that  Ja^  =  J3*q  +  ja^ 
(see  (2.4))  and  therefore  -^N  =  -^Nj0  +  ^,1’  where  ^.0  =  PN  and 

■£*N  ^  =  PNJa*jPN.  q  is  a  bounded  linear  operator  on  Z  (see  [3],  Proof 

of  Theorem  3.1).  We  only  need  to  estimate 

|-Q4,lZl“a1si,lz2|z  -  |J^lPNZl“'aiPNZ2|z 

|  f  (4>N) -f  (^N)  |  i  L(^o|<|>N(-r;j)-i|)N(-rj)  [  +  |4>n-^n|2> 

—  L(m-r<s<0  l<|,N(a)-*N(s)|  +  |PNZl-PNZ2|z)* 
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where  z^z,,  €  Z  and  P^  =  (<i>N(0)  ,<J>N)  ,  PNz2  =  (^(0),^).  Since  ZN  is 

finite  dimensional  all  norms  on  Z„  are  equivalent.  Therefore  there  exists 

N 

a  constant  0  >  0  such  that 

N 

!"^N,lZl  "  J3^N,1Z2|z  -  °n|Z1"  Z2  |  Z 

In  order  to  prove  part  c)  we  recall  that  dissipativeness  of  -  toi 

implies  the  existence  of  (I  -  ^  on  i^(I  ~  for  ^  ^  (0,^). 

Assume  that  ^(1  -  A.,0/^)  ■  Z.  We  take  y  f  and  put  z  =  (1  -  ^y* 

Then  z  =  z^  +  z2,  where  z^  €  and  z2  €  Z^,  and 

zx  +  z2  -  ^-Q^j(z^+Z2)  =  y  €  Zj^ 

implies  z2  =  0,  i.e.,  z  £  ZN>  Note,  that  ja^Z  ^  ZN> 

Proposition  3.1. 

a)  _Q^n  generates  a  semigroup  TN(»)  of  type  w  on  z  (i.e., 
properties  (i)  -  (iv)  of  Lemma  2.1,  b)  hold  for  TN(’)  with  M  =  1) , 

T  (t)z  =  lim  (I  -  ^0"nz,  z  £  Z, 
n-*» 

uniformly  with  respect  to  t  in  bounded  intervals. 

b)  TN(t)ZN  c  ZN  for  t  >  0,  H-1,2,...  . 

c)  For  any  z  £  Z  the  function  u(t)  =  T^(t)z,  t  _>  0  is  continuously 
differentiable  for  t  _>  0  and 

u(t)  -  -Q^u(t),  t  _>  0, 


u(0)  ■  z. 


'3.1) 
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Proof.  Since  by  Lemma  3.1,  b)  S3^  is  globally  Upschitzean  on  Z, 
equation  (I  -  Xj^)  z  =  y  has  for  any  y  6  Z  a  solution  z  £  Z  provided 
X  is  sufficiently  small,  i.e.,  <&(I  -  XiaO  =  Z  for  those  X.  This 

together  with  dissipativeness  of  -  ml  (Lerana  3.1,  a))  shows  that  we  may 

use  the  Crandall-Liggett  Theorem  ([8;  Theorem  I])  in  order  to  get  a).  Part  b) 
then  is  an  immediate  consequence  of  Lemma  3.1,  c) .  It  is  clear  that  equation 
(3.1)  for  any  z  €  Z  has  a  unique  solution  on  [0,°°)  which  is  continuously 
differentiable.  By  Theorem  II  of  [8]  we  get  u(t)  =  T^(t)z. 

In  the  proof  of  Proposition  3.1,  a)  we  have  seen  that  for  each  N 
iSd  -  X_£*^)  =  Z  for  X  sufficiently  small.  For  our  convergence  proof  we 
shall  need  that  this  range  condition  holds  uniformly  with  respect  to  N. 

Proposition  3.2.  Let  to  be  the  constant  of  Lemma  3.1,  a) .  Then 

^(1  -  Xj/N)  *=  Z 

for  all  X  6  (0,—)  and  all  N. 

ur 

Proof.  In  the  proof  of  Proposition  3.1,  a)  we  have  seen  that  for  any  N 

there  exists  a  X^  >  0  such  that  -  XjQ*')  =  Z  for  all  X  €  (0,XN). 

Then  for  X  <  min(X„ )  we  have 
—  N  to 

z  =  mi  -  xj^N)  =  ^(d-Xto)d  -  -  wI>>> 

'  mx  -  dh'^N  -  • 

i.e.,  ^(1  -  y(_Q^  -  tol))  ■  Z  for  y  sufficiently  small.  Since  -  wl 
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is  dissipative,  we  immediately  (see  for  instance  {6;  p.73])  get 

j^(I  -  -  oil))  *  Z  for  all  U  >  C  or  equivalently  &(I  -  X_Q^)  =  Z 

for  all  X  £  (0,i). 

’id 

By  Proposition  3.1  we  obtain  a  sequence  of  ordinary  differential 
equations  on  the  finite  dimensional  spaces  Z^, 

zN(t)  -j/NzN(t),  t>0,  (3.2) 

where  denotes  the  restriction  of  to  Z„.  For  any  z  €  Z  the 

N  N  N 

solution  of  (3.2)  with  initial  condition  CO)  =  PNz  is  given  by 

zN(t)  =  T^(t)PjjZ.  In  order  to  prove  that  zN(t)  T(t)z  as  N  -*  °°  we  shall 

use  the  following  nonlinear  version  of  the  Trotter-Kato  Theorem: 

Theorem  3.1.  Let  N  »  1,2,...,  and  be  single-valued  operators 

on  a  Banach  space  X  such  that  domJ#N  o  dom  for  all  N  and  domi#  =  X. 
Moreover,  assume  that  the  following  conditions  are  satisfied: 

(i)  There  exists  a  Xq  >  0  such  that 
-  \0)  -  -  X  4^)  m  X 

for  N  =  1,2,...  and  all  X  €  (0,Xq). 

(ii)  There  exist  real  constants  (D^,  N  *  1,2,...,  and  id 
such  that  the  sequence  {id^}  is  bounded  above  and 
-  wNI  and  -  <Dl 
are  dissipative  in  X. 

(iii)  There  exists  a  subset  &  of  dom  such  that 
(I  -  X0)S*  ■  X  for  X  sufficiently  small  and 
^Nx  as  N 

for  all  x  6 


GO 


Let  and  T(*)  denote  the  semigroups  generated  by  and 

respectively.  Then  for  all  x  €  X 

lira  T„(t)x  =  T(t)x 
N-*» 

uniformly  with  respect  to  t  in  bounded  intervals. 

This  theorem  can  quite  easily  be  extracted  from  [7].  Using  the  same 
conclusions  as  in  the  proof  of  Theorem  4.1  in  [7]  one  first  shows  that  there 
exists  a  X  >  0  such  that  (I  -  (I  -  Xj#)^  as  N  -►  <*>  for  all 

x  €  X  and  all  X  €  (0,X).  The  rest  of  the  proof  is  the  same  as  the  proof 
of  Theorem  3.1  in  [7]. 

Theorem  3.2,  Let  be  an  approximation  scheme  for  equation  (2.1) 

and  assume : 

(i)  P^z  -*■  z  for  all  z  €  Z. 

(ii)  There  exists  a  k  >  1  such  that  <j>  £  C  implies 
4>n(0)  -*•  4>(0)  and  |<f>N  -  ^ 1 2  *  0 
as  N  <j>N  €  is  defined  by 

PN(<t>(0),<J>)  -  (<|>N«)),<J>N). 

Then  for  any  z  €  Z 

lira  TN(t)PNz  -  T(t)z 
N-x» 


uniformly  on  bounded  t-intervals. 
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Proof.  By  assumption  (i)  it  is  sufficient  to  prove  lim  TN(t)z  =  T(t)z. 

N-KD  N 

Conditions  (i)  and  (ii)  of  Theorem  3.1  are  satisfied  because  of  Lemma  3.1,  a) 

and  Proposition  3.2.  The  verification  of  condition  (iii)  in  Theorem  3.1  is 

done  in  the  same  way  as  in  the  linear  case  choosing  =  .5^  (see  [5], 

Proof  of  Theorem  3.1).  Note,  that  under  condition  (ii)  we  have 

sup  |<J>  (s)  — (s)  |  ->  0  and  therefore  by  (H2)  f(<J>  )  -*■  f(<J>)  as  N  ■*  ». 

-r<s<0  N  N 


Remark.  The  proof  of  Theorem  3.1  in  [5]  uses  only  Q)  =  and  not 

g,  m 


h" 


a1;  = 


Let  =  dim  ZN.  We  choose  a  basis 
(8^(0)  ,3j),3j€  and  define  the  n  x  k^  matrix 


3N 

*1' 


(3i » •  •  • 


*>■ 


rn 


of  Z 


N» 


where 


If  we  put  gN  =  (3N(0),gN),  then  for  any  z^  €  ZN  we  have 

ZN  =  &N“n  =  (&N(0)c*N.BNaN), 

^N 

where  6 1R  is  the  coordinate  vector  of  z^.  The  same  calculations  as 
in  [5]  show  that  for  (ti,4>)  €  Z 

PN(n,<t>)  -  8%*, 


where 


is  the  solution  of 


■  kjj  (r) ,  4> ) . 
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The  *  k^-matrix  QN  and  the  kjj-vector  h  (n.4>)  are  given  by 


Qn  =■  6N(0)T6N(0)  + 


(s)g(s)ds, 


[°  8N<9)TeNi 
J  -r 

h^M)  -  SN(0)Tn  +  |°  eN(s)T4>(s)g(s)ds. 


The  matrix  representation  ^  of  q  -  PN**0PN  restricte<l  to  ZN  is 

given  by 

N.o  ”  qn  hn* 

where  =•  |  BN(s)T  |3N(s)g(s)ds. 

For  zN  €  ZN  with  coordinate  vector  we  denote  the  coordinate 

vector  of  .Q^  ^  *  Pjj.Qf'P^  with  ^^jj)  •  From  izN  ”  PN^f^NcV’0* 

we  get  that  FN(aN)  is  the  solution  of  *  hN(  (f  (S**^)  ,0) )  - 

6N(0)Tf(3NaN),  i.e., 

FN(0lN)  " 


zN(t)  is  a  solution  of  (3.2)  with  zN(0)  “  pN(n,<}>)  if  and  only  if  the 

'■'N 

coordinate  vector  w^(t)  of  z^(t),  zN(t)  ■  3  w^ft),  is  a  solution  of 

WN^  "  AN,OwN^t^  +  FN*wN^t^* 

WN<0)  -  djj,  where  3^  -  PN(n,$). 


(3.3) 


15. 


By  Theorem  3.2  we  have 


;N 


lim  8  «N(t) 
N-h»  N 


(x(t;n,4>),  xt c ) 


(3.4) 


uniformly  on  bounded  t-intervals. 

4 .  Spline  Approximation 

A  realization  of  the  scheme  presented  in  Section  3  can  be  obtained  by 
using  subspaces  of  spline  functions.  In  order  to  do  so  take  a  sequence  of 
partitions  of  [-r,0]  with  meshpoints  t^,  N  ■  1,2,...,  such  that 


max  1 1 .  -  t  .  |  0  as  N  -*  “>  and  max 

j  J  J 


jX  "  tj-ll/m|nltj  "  t ”_i  I  i  6  <  00 . 
N  -  1,2,...,  for  some  constant  8  >  0.  For  instance  we  may  take 


N  N 


ej  •  -jf.  j  *  o, . . . ,n. 


or 


„N  4  rk+l“rk  .  -  1  4  ,  V 

^kN+j  "  ”  j  jg  »  ^  m  0, . . .  ,m  •  1,  j  "  1 , . . .  ,N. 


In  the  second  case  -r^t  h  *  0,...,m,  is  always  a  meshpoint  of  the  partition. 

Theorem  4.1.  Let  Z„  be  the  subspace  of  all  elements  (<$>(0),$)  in  Z  such 

- ■  N 

that  $  is  a  first  order,  a  cubic  Hermite  or  a  cubic  spline  with  knots  at 

tj,  respectively,  and  choose  the  corresponding  approximation  scheme  {Z^.P^.j*^} 

for  equation  (2.1).  Then  for  any  z  €  Z, 


lim  TM(t)PNz  -  T(t)z 
N-~>  N  N 
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uniformly  on  bounded  t-intervals,  where  z..(t)  -  T.,(t)P.,z  is  the  unique 

N  N  N 

solution  of  (3.2)  with  z„(0)  •  P„z. 

Proof.  We  have  to  verify  conditions  (i)  and  (ii)  of  Theorem  3.2.  This  is 
done  in  the  same  way  as  in  the  linear  case  using  the  nice  convergence  properties 
of  spline  functions  (see  [5],  Proof  of  Theorem  4.1,  or  [9],  Section  4). 

For  calculations  one  chooses  B-splines  as  basis  elements  of  (see 

for  instance  (14],  Sections  2.1,  3.1  and  4.1).  The  matrices  Q  ,1L^  are 
calculated  as  indicated  at  the  end  of  Section  3.  Of  course,  system  (3.3)  is 
always  solved  in  the  form 

Qnwn(0  “  +  8N(0)Tf  (3NwN(t))  . 

The  cubic  spline  scheme  was  used  for  the  following  initial  value 

x(t)  -  -x(t-l) (4-x2(t)) ,  t  ]^0,  (4.1) 

x(t)  •  cos  j  t,  -1  <  t  <  0.  (4.2) 

Nussbaum  proved  [13;  Corollary  2.1]  that  equation  (4.1)  has  a  unique  slowly 
oscillating  solution  Xg(t)  such  that  Xg(t)  >  0  on  (0,21,  Xg(-t)  ■  -xQ(t) 
and  Xg(t+2)  ■  -Xg(t)  for  all  t.  Moreover,  { Xp(t ) j  <  2  for  all  t.  The 
solution  of  (4.1),  (4.2)  tends  to  Xg(t+c)  for  some  c  €  IR. 

In  Table  1  we  give  the  values  for  x(t)  on  the  interval  (0,5], 


17. 

which  were  obtained  by  stepwise  integration  of  (4.1),  (4.2)  using  a  standard 

U 

fourth  order  Runge-Kutta  scheme,  and  the  values  of  B  (0)w  (t),  N  -  16, 
where  wN(t)  is  the  solution  of  (3.3)  for  n  *  1,  4(t)  -  cos  j  t.  A  is  the 
difference  x(t)  -  B^^(0)w^^(t) . 
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TABLE  I 


t 

x(t) 

B16(0)v16(t) 

A.10-5 

0.0 

1. 

1. 

0 

0.2 

0.90365 

0.90365 

0 

0.4 

0.59385 

0.59384 

1 

0.6 

0.04893 

0.04890 

3 

0.8 

-0.63787 

-0.63791 

4 

1.0 

-1.23865 

-1.23869 

4 

1.2 

-1 . 60894 

-1.60897 

3 

1.4 

-1.77853 

-1.77855 

2 

1.6 

-1.82899 

-1.82901 

2 

1.8 

-1.78650 

-1.78653 

3 

2.0 

-1.56826 

-1.56829 

3 

2.2 

-0.89005 

-0.89006 

1 

2.4 

0.40275 

0.40279 

-4 

2.6 

1.45999 

1.46003 

-4 

2.8 

1.85918 

1.85920 

-2 

3.0 

1.96289 

1.96289 

0 

3.2 

1.98664 

1.98664 

0 

3.4 

1.98925 

1.98926 

-1 

3.6 

1.97633 

1.97634 

-1 

3.8 

1.90934 

1.90935 

-1 

4.0 

1.61011 

1.61013 

-2 

4.2 

0.62214 

0.62218 

-4 

4.4 

-0.88263 

-0.88258 

-5 

4.6 

-1.70637 

-1.70636 

-1 

4.8 

-1.93456 

-1.93456 

0 

5.0 

-1.98423 

-1.98422 

-1 

! 

i 
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